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^ ■ ABSTRACT 

f~i , Spectroscopic modeling of Type II supernovae (SNe) generally assumes steady-state. Fol- 

Qh' lowing the recent suggestion of Utrobin & Chugai, but using the ID non-LTE line-blanketed 

Q \ model atmosphere code CMFGEN, we investigate the effects of including time-dependent 

5^ . terms, which are generally neglected, that appear in the statistical and radiative equilibrium 

' equations. We base our discussion on the ejecta properties and the spectroscopic signatures 

^ \ obtained from time-dependent simulations, investigating different ejecta configurations (slow, 

' i standard, and fast), and covering their evolution from one day to six weeks after shock break- 

^ ' out. Compared to equivalent steady-state models, our time-dependent models produce SN 

, ejecta that are systematically over-ionized, affecting helium at one week after explosion, but 

■ ultimately affecting all ions after a few weeks. While the continuum remains essentially un- 

QQ I changed, time-dependence effects on observed spectral lines are large. At the recombination 

. epoch. Hi lines and Nal D are considerably stronger and broader than in equivalent steady- 

' state models, while Call 8500A is weakened. If time dependence is allowed for, the Hel lines 

\ at 5875A and 10830A appear ^3 times stronger at one week, and Hel lOSBOA persists as 

I i a blue-shifted absorption feature even at 6 weeks after explosion. Time dependence oper- 

' ates through the energy gain from changes in ionization and excitation, and, perhaps more 

\ universally across SN types, from the competition between recombination and expansion, 

' which in-turn, can be affected by optical-depth effects. Our time-dependent models compare 

well with observations of the low-luminosity low-velocity SN 1999br and the more standard 
SN 1999em, reproducing the Ha line strength at the recombination epoch, and without the 

^ ■ need for setting unphysical requirements on the magnitude of nickel mixing. 
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I INTRODUCTION 

Because of radiative cooling and the fast expansion of the explod- 
ing mantle of the progenitor massive star, photospheric-phase Type 

II SN spectra evolve rapidly. At shock breakout the spectral en- 
ergy distribution (SED) peaks in the far-UV and X-rays. Subse- 
quently the SED centers in the UV, then in the optical, and later 
in the infrared until the object fades from view. From a radiative 
transfer perspective, the cooling induces a recombination to lower 
ionization stages that impose a strong blanketing effect on the en- 
ergy distribution. A few days after explosion, the ejecta radiate a 
UV-dominated SED, with Hell 4686A (only observed at a few days 
after shock breakout; Dessart et al. 2007), Balmer/Paschen lines 
of Hydrogen, Hel5875A, Hel 10830A, some multiplets of Oil and 
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Nil multiplets at 4600A (see Dessart & Hillier 2006a; Baron et 
al. 2007), Nil 5400A, and a few isolated resonance lines such as 
MgIl2802A and Aim 1859A. As the ejecta continue to cool, hy- 
drogen eventually recombines, with a contemporaneous enhance- 
ment in line-blanketing due to the switch from Felll to Fell. In 
Dessart & Hillier (2006) and Dessart et al. (2007), we success- 
fully modeled these epochs for SNe 1999em, 2005cs, and 2006bp, 
using steady-state non-LTE CMFGEN models and a power-law 
density distribution with an exponent of ten (the immediate post- 
breakout phase requires higher values for this exponent). These 
spectroscopic analyses deliberately focused on the first ~45 days 
after shock breakout, since, beyond that time, we encountered se- 
vere difficulties in reproducing the hydrogen lines. Specifically, at 
late times, synthetic line profiles were sizably narrower than ob- 
served, suggesting that line formation regions predicted by CMF- 
GEN were confined to velocities/radii that were too small. Despite 
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a continued success with reproducing the overall continuum energy 
distribution, CMFGEN failed to reproduce important line profiles, 
whenever hydrogen started recombining at and above the photo- 
sphere. 

The Ha problem has not been clearly emphasized in the 
literature, where one can see a great disparity in model atmo- 
sphere assumptions and agreement between theoretical predictions 
and observations, but no direct link to physical/numerical issues. 
In the eighties and early nineties, the recognized importance of 
non-LTE effects confronted the strong limitations of computer 
technology, so that only a few species were treated in non-LTE 
(i.e., usually hydrogen and helium), while the metals responsi- 
ble fo r line blanketing were treated in LTE. lEastman & KirshneJ 
( Il989i) followed this approach to model the first ten days of SN 
1987 A, and did not encoun ter obvious difficulties with hydrogen 
lines. ISchmutz et al. using an approximate non-LTE tech- 

nique, had, on the contrary, great difficulty reproducing any of the 
Balmer lines, suggesting clumping as the culprit. iHoflichl l ll988l) 
reproduced the SN 1987 A spectral evolution and the hydrogen 
lines, over many months, using a large " turbulent" velocity. In 
the more sophisticated non-LTE CMFGEN teillier c& Miller" 1998; 
jDessart & Hillier 2005a) models presented in Dessart & Hillier 
l l2006ah . we found that decreasing the turbulent velocity weakened 
line-blanketing effects and increased, although only modestly, the 
strength of hydrogen lines. However, beyond 40 days after explo- 
sion, this tuning had no longer any important influence on the hy- 
drogen lines. 

The non-LTE model atmosphere code PHOENIX pre- 
dicts strong B almer line s at the hydro gen recombination epoch 
jMitchelletall 2001; Ba ron et al.l [20031) , but this may stem from 
their adoption of non-thermal ionization/excitation due to ^^Ni at 
the photosphere, sometimes just a few days after ex plosion and in 
mass shells moving at > 10000 km s"^in S N 1987 A jMitchell et al] 
I2OOI I) or ~9000km s^^in SN 1993W ( iBaron et al.l l2003h . Hy- 



drodynamical simulations of core-collapse SNe predict that ^^Ni 
has velocities of at most ~4000km s~^, th e nickel fingers be- 
ing strongly decelera t ed at the H/He interface jFrvxell et al.lll99ll ; 
iKifonidis et alll2000l l2003h . The magnitude of this disagreement 
extends far beyond the uncertainties of explosion models and sug- 
gests a genuine incompatibility. 

Interestingly, Hi lines remain strong for months in all Type 
II SNe, in objects as diverse as the "peculiar" SN 1987A, the 
"plateau" SN 1999em, and the "low-luminosity" SN 1999br. CMF- 
GEN models computed for these objects were unable to reproduce 
Balmer lines after ~4 days in SN 1987 A, ~40 days in SN 1999em, 
and ~20 days in SN 1999br, all coincident with hydrogen recom- 
bination in the ejecta. These three SNe have very different inferred 
ejecta properties and observed light curves. SN 199 9br even synthe- 
sized an order of magnitude less ^^Ni than average dPastorello et al.l 
|2004) for the Type II class. The only common property between 
these objects, which is connected to the Ha problem, is the recom- 
bination of the ejecta to a lower ionization state at the correspond- 
ing epoch. 

Recently, lUtrobin & Chugail ( l2005l UC05) proposed that the 
effect of time-dependence, and the energy associated with changes 
in ionization/excitation, lead to a strong Ha line profile in SN 
1987 A during the recombination epoch. They also found that bar- 
ium lines were affected, and that, with their more consistent ap- 
proach, BaIl6142A could be fit ted using the LMC metallicity 
value. The steady-state models of iMazzali et al.l ( Il992h supported 
instead an abundance enhancement of five. 

Time dependence has been invoked in the past by 



iFransson & Kozm al l ll993h to explain the late-time ligh t curve of 
SN 1 987 A, and the theoretical study of lPinto & Eastmatj JioOO^Ibl) 
showed that time dependence in the radiation fiel d had a criti- 
cal imp act on the radiative transfer in Type las. Pi nto & EastmanI 
( l2000al lbl). however, treat the material in LTE, i.e., do not solve 
the rate equations, focusing instead on the time-dependent dijfu- 
sion o f photons through a n optically thick Type la SN ejecta. Sim- 
ilarly, 'Kas en et"al] 1I2OO6') neglect explicit time dependence in the 
rate equations, computing the ionization and excitation state of the 
medium in LTE. Thus, although there is at present growing interest 
in accounting for time dependence in the radiation field, the often- 
used expedient of LTE, to maintain low CPU costs, has forced the 
neglect of both non-LTE and time-dependence in the level popu- 
lations. On e exception to t hese time-dependent LTE approaches is 
the work o f lHoflichl ( l2003l) . who treats time-dependence in the rate 
equations but, to our knowledge, has not discussed the associated 
effects on Type II SN spectra. In the present study, we investigate 
thoroughly the effects of time-dependence in the rate equations and 
their impact on inferred ejecta properties. Note that time depen- 
dence in the radiation field is accounted for by adjusting the base 
luminosity so that the emergent synthetic flux matches the observed 
flux, using the bolometric -light evolution of SN 1999em as a guide 
jPessart & Hilliej2006ah . 

Here, we report the salient features of several time-dependent 
non-LTE CMFGEN simulations, covering the early evolution for a 
range of Type II SN ejecta. We confirm the results of UC05 that 
time-dependence induces an over-ionization of the recombining 
ejecta and that it solves the Ha problem. However unlike UC05, we 
self-consistently solve the radiation transfer equation — the cou- 
pling between the level populations and the radiation field is calcu- 
lated and fully allowed for. Our more comprehensive study of the 
full spectrum further predicts that all lines, not just those of hydro- 
gen, are significantly affected. The ionization of the ejecta and its 
evolution are so strongly modified that we predict even He: 10830A 
many weeks after explosion. 

In the next section, we present out treatment of time- 
dependence in CMFGEN, which focuses here on the terms ap- 
pearing in the statistical and radiative equilibrium equations (an ap- 
pendix also provides further details). In §2.5, we present the various 
model calculations performed to illustrate our discussion. In §3, we 
present our results, discussing the effects of time dependence on the 
ejecta properties as well as the associated spectroscopic signatures. 
We then present in §4 a comparison with a few representative obser- 
vations, focusing on the well observed Type II Plateau SN 1999em 
and the low luminosity Type II SN 1999br. In §5, we discuss the 
implications of such time-dependent effects on our understanding 
and on our modeling of photospheric-phase Type II SN spectra, be- 
fore giving our conclusions in §6. Note that the salient features and 
key results presented here are also available in a concise, "letter", 
format in .Dessart & Hillier, ( .2006h) . to which we refer the hurried 
reader. 



2 TREATMENT OF TIME DEPENDENCE 
2.1 Statistical Equilibrium Equations 

The general form of the time dependent statistical equilibrium 
equation is 



+ V.(niti) = {njRji - UiRij) 



3 



(1) 
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where v is the velocity, Ui is the population of state i, and Rij is 
the transition rate from state i to state j. The forms of Rij and Rji 
are standard, and can be found, for example, in Mihalas (1978). 
For simplicity we have not distinguished between the different pro- 
cesses (collisional, bound-free, and bound-bound) entering these 
rates. 

In massive stars, the winds can be approximated as stationary. 
In such cases we consider a Eulerian formulation, and take dn/dt 
to be zero. The V .{jiiv) term is the advection term, and is usually 
neglected, however it can become important in the outer regions of 
the wind, since the ratio of the recombination time, l/an^, to the 
flow time, r/v, scales as the radius, r (here, a is the recombination 
coefficient). 

Supernovae are obviously time dependent, and in this case it 
makes sense to use a Lagrangean formalism. Using the continuity 
equation 



^+P(V.«) = 0, 

the statistical equilibrium equations become 

Dui/p 



Dt 



where D/Dt is the comoving derivative. For a Hubble flow, 

1 Dr'^m 



r3 Dt 



E 



{rijRji — TiiRij) 



(2) 



(3) 



(4) 



which is the form we used in these calculations. 

To solve these equations we use implicit first order differenc- 
ing in the Lagrangean frame. Thus the equations to be solved be- 
come 



(ni)k — (wi)fc_i(rfc-i/rfc)^ 
At 



j 



Tli Rij ) (, 



(5) 



where k — l,k refer to values at consecutive time steps but with the 
same comoving coordinate. For a Hubble flow, we can simply use 
the velocity as the Lagrangean co-ordinate. The solution of these 
equations proceeds in an identical fashion to that used to solve the 
steady-state statistical equilibrium equations (Hillier 1987; Hillier 
1990; Hillier & Miller 1998); the only distinction is that we have 
a "source" term that comes from the populations at an earlier time 
step. 



2.2 Energy balance 

The energy equation is 



De PDp ^ ^ . 



(6) 



where e is the internal energy/unit mass, the opacity, r;^ the 
emissivity, and J,j the mean intensity, and we have ignored other 
forms of energy deposition (such as nuclear decay), e can be written 
in the form 



e = ck + ei , 
where 

3A;r(n + Ue 



and 



2pmn 



(7) 



(8) 



6/ 



EriiEi 
pmn 



(9) 



In the above, n is the total particle density (excluding electrons). 
Tie is the electron density, and Ei is the total energy (excitation 
and ionization) of state i. This equation is solved implicitly, via 
linearization, and is no more difficult to treat than the regular con- 
straint of radiative equilibrium. 



2.3 Are the time-dependent terms important? 

Below we examine when the time dependent terms in the statistical 
equilibrium equations will play an important role in determining 
the radiative transfer in SN envelopes. To do this, we need to com- 
pare the time it takes to replenish the population of individual levels 
with the flow time. As noted previously, the flow-time scale is sim- 
ply r/v, and for a Hubble flow this is texp, where texp is the time 
since the explosion. 

We will only examine the ionization equilibrium in detail, 
since it is the ionization equilibrium that will be most affected by 
the time-dependent terms - statistical equilibrium of excited levels 
will generally be satisfied because the radiatives rates tend to be 
large {A > 10® s~^), and the transitions have low (or moderate) 
optical depths. 

The simplest recombination timescale can be found by esti- 
mating the time it takes, at fixed electron density, for all ions to 
recombine. Thus the recombination timescale for hydrogen, tree is 
simply given by 



tree 



x(H+) 



n(H+)nea 



Using a = as - 

dOsterbrockll 19741) . gives 

4.47 X 10'' , 
^rec = days. 



(10) 

2.59 X 10~^^cm^s"^ at 10,000K 



(11) 



We use as since in SN conditions every direct recombination to 
the ground state will generally be followed by an ionization. The 
recombination time scale for Hel is very similar to that of Hi. Since 
the comoving density is oc , the ratio of recombination time to 
the flow time, at a given velocity, scales as . Thus the ionization 
must eventually become frozen, and first at high velocities. 

The previous discussion refers to the case when we can ig- 
nore optical depth effects and the existence of metastable states. 
When optical depth effects are important, and/or there exist low ly- 
ing metastable states, it is more difficult to ascertain the relevant 
recombination time scale. In such cases, it is the net flow of elec- 
trons to the ground state that will set the ionization equilibrium. 

As an illustration consider a two level atom with continuum. 
The relevant equations are 

In 

and 

Dn2 
Dt 
and 



= n2-R2i - niRi2 , 



n^nia — n2P2i + niRi2 ~ n2R2 



(12) 



(13) 



n-i + n2 + ni = n 



(14) 



Here nj is used to denote the ion population, P21 is the photoion- 
ization rate, R21 refers to the combined rate for processes which 
couple level 2 to 1 (bound-bound, collisional, two photon decay). 
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while Ri2 refers to processes that excite an electron from level 1 
to 2. We assume that the ground state continuum is in detailed bal- 
ance. For a strict two-level atom a — ai2 but for these discussions 
we can take a — ub- 

In the absence of the time-dependent terms, the ground state 
equation reduces to 



n2-R21 ~ n\R\2 



(15) 



In the case of Hel the main decay route from the S state is two 
photon decay (with A = 52 s~^), while coUisional excitation might 
be the main reverse mechanism. 

When we include the time-dependent terms, the equation for 
level 1 can reduce to 



Dt 



— n2-R21 : 



and that for level 2 becomes 



'^2-P2/ = rieTija , 



(16) 



(17) 



since the decay term, for a metastable level or for a level decaying 
through an optically thick transition, can be small. The effective 
recombination time scale, t^ff, is (since Dni/Dt — —DX^/Dt) 
thus 



'^eff ^ 



Dni/Dt 

12-R21 
^ 

irec^. 



(18) 
(19) 
(20) 

(21) 



Thus, if the photoionization rate P21 is much greater than the 
decay rate R21, the effective recombination time scale will be con- 
siderably longer than the classic recombination time scale. 

The expression for P21 is 



dv , 



(22) 



where o is the photoionization cross-section and h is Planck's 
constant. Assuming Ji, = W Bi, (a diluted blackbody), a = 

f7"o(!^o/i^)^, and /ii^o/fcr >> 1, we have 



-21= 5.8xlO«W^(^)(^) 

1 



exp{—hVo/kT) s 



For Hel 2s ^S, we have yl(two phot) = 52s"\ Uo = 0.960 x 
lO^'^Hz, ao = 8.90 x 10~^*cm^, and, assuming W = 1, and 
T = lO* K, we find 

tgff = 9.2 X 10^ tree • 

Thus, it is apparent that the effective recombination time scale is 
much longer than the classic recombination time scale, and could 
easily be comparable to the flow time scale. 

We appreciate that the true situation is more complex than 
that discussed above. For Hel the 2s and 2p^S state are both 
metastable, and there is a collisonal coupling between the sin- 
glet and triplet terms. The 2s state can decay directly to the 
ground state, for example, by two-photon decay, by collisional de- 
excitation, and by collisional excitation to the 2p state followed 
by a radiative decay. Further, the effective recombination rate will 
substantially exceed the classic recombination time scale only if 



the resonance line is optically thick. Despite these caveats, it is ap- 
parent that the time-dependent terms might be important even at 
the photosphere, and must be included. A simple comparison be- 
tween the flow time scale and the classic recombination time scale 
can greatly underestimate the importance of these time-dependent 
terms. 

In the Appendix we examine, in detail, the various processes 
that influence the ionization of helium at a representative depth, 
and in a typical model, and where the flow time is significantly 
greater than the classic recombination time. In that appendix we 
show that the combination of optical depth effects, and the exis- 
tence of metastable levels, leads the time-dependent terms to cause 
a delay in the recombination of He+ to neutral He. 

For Hi, the 2s state is metastable, but in atmospheric calcu- 
lations (as is done here) it is usually assumed to be in LTE with 
respect to the 2p state (i.e., the 2s and 2p states have the same de- 
parture coefficients). Thus the relevant time scale is set by radiative 
decay from level 2p to Is, and, hence, by the Lya optical depth, as 
was previously noted by UC05. Other processes or configurations 
(e.g., overlapping lines, continuum opacity) can also be important 
since they can destroy a Lya photon, facilitating a "real" recombi- 
nation to the ground state. 

Finally we stress, as apparent from the above discussion, that 
the time-dependent term does not directly change the ionization 
equilibrium (unlike the situation at low electron densities). Rather, 
it prevents the decay of electrons from the n = 2 states (for H I and 
He I) to the ground state. The ionization of H I and He I can then be 
maintained by ionizations from these excited states. 

2.4 Radiative Transfer 

Explicit in our calculations is the assumption that we can solve 
for the level populations and temperature structure neglecting the 
relativistic and the time dependent terms in the radiative transfer 
equation. While it is moderately easy to include the relativistic 
terms (e.g., Hauschildt 1992) the time depende nt term is equally 
import ant and must be included for consistency jPinto & EastmanI 
l2000allbl) . In hydrodynamic models, both terms are routinely (and 
must be) included, both for the hydrodynamics and computing the 
light curve. 

However for spectral calculations it is usual to follow a stel- 
lar approach and assume that we can compute a static atmosphere 
by specifying the luminosity via the diffusion approximation at the 
lower boundary. Studies of the importance of the relativistic terms 
have been done, for example, by Hauschildt et al. (1991), or by Jef- 
fery (1993), and shown to be important. What is not clear, however, 
is what effect the terms will have on the distances based on the Ex- 
panding Photosphere Method (Dessart & Hillier 2005b), something 
that we are currently investigating. Very preliminary work suggests 
that the inclusion of relativistic terms, while producing spectra that 
are very similar, could lead to changes in EPM distance by up to 
10%. However, we do not expect that the inclusion of these terms 
will effect the conclusions in this paper; the photospheric temper- 
atures are similar and the importance of the time-dependent terms 
arises primarily through atomic physics, and should not be sensitive 
to radiative transfer affects. 

2.5 Presentation of models 

We adopt a different procedure to mode l SN sp ec tra from that em- 
ployed previously in iDessart & Hillied i l2005a llbl l2006ah by mod- 
eling the time evolution of a larger portion of the ejecta, and in a 
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Table 1. CMFGEN ejecta characteristics. From left to right, the column refers to the model name, the initial base radius Rq, the base velocity vq, the initial base 
density po, the total mass of the ejecta covered by the CMFGEN grid .A/(,jccta. the initial and final times ti and tj, the initial and final base electron-scattering 
optical depths t^sA ™d r^s f, assuming full and uniform ionization (kcs = 0.30), the initial and final base electron-scattering optical depths t^s i and t^^ f 
computed with CMFGEN including time-dependent terms in the statistical and radiative equilibrium equations, and, in the last column, the final base electron- 
scattering optical depth computed with CMFGEN under the assumption of steady-state. Note how the value of t^s f diminishes from the assumption of fixed and 
full ionization, CMFGEN with time-dependence, and CMFGEN without time-dependence. 



Ro fO PO Afcjccta *i *f T"cs,i Tcs,f fos,i T"cs,f f"cs,f 
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D/DT no D/DT 
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Figure 1. Synthetic light curves for the baseline model (A), the slow (B) and the fast (C) version of model A, and the "early" model (D), which covers the first 
week after explosion. These models correspond to standard (models A and C) and underluminous (models B and D) Type II SNe, and reproduce approximately 
the evolution of their bolometric luminosity and color, with a bolometric luminosity that decreases and an SED that peaks further and further to the red with 
time. Note that in some models, the emergent luminosity shows a small rise at the end of the sequences, resulting from a slight inadequacy in the imposed flux 
at the inner boundary. [See the electronic edition of the Journal for a color version of this figure.] 



Lagrangean spirit. Each model is stepped in time using the implicit 
method outlined above. We have tested the dependence of our re- 
sults on the adopted time increment. Two sequences evolved with a 
1 day increment or a 2 day increment, but starting and ending with 
the same model input parameters, yield an SED at the final time 
that agree to better than 0. 1% on average, departures peaking at 1% 



at certain UV wavelengths (coinciding with strong fines but asso- 
ciated with spectral regions of negligible flux). In practice, we in- 
crease the time increment (sfightly) with the time elapsed since ex- 
plosion, using a fraction of a day within the first week, but ~3 days 
at a few weeks after shock breakout. For simplicity, we choose the 
same minimum and maximum velocities for all models in a given 
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Figure 2. Left: Radial variation of the ionization fraction of Hydrogen (Hll: solid line; Hi: Broken line) for the baseline model A with time dependence (black) 
and without (red), and at 48.7 days after explosion, but otherwise having identical parameters. Note the over-ionization of hydrogen at large radii, not present 
in its neutral state despite the very late time, akin to a frozen-in ionization. Symbols give the positions of the CMFGEN adaptive-grid points, for each model 
but for only one curve. Note the presence of the well-resolved recombination front near the base, deeper-in for the steady-state CMFGEN model. Right: Same 
as left, but for the ionization fraction of helium (Hell: solid line; Hel: Broken line). Notice the same frozen-in ionization at large radii. [See the electronic 
edition of the Journal for a color version of this figure, and the text for discussion.] 



sequence. Because we neglect the relativistic and time dependent 
terms from the transfer equation, the base luminosity is a model 
parameter only, and does not reflect the true physical base lumi- 
nosity. Implicit in our method is the assumption that the observed 
spectrum primarily depends on the local photospheric conditions 
as set by the observed luminosity. We use the same approach in 
the steady-state modeling, although in that case we only need to 
model the atmosphere to a depth where we can utilize the diffusion 
approximation, usually corresponding to a Rosseland optical depth 
of a few tens. Here, the initial inner radius of the computation is 
chosen so that it remains optically thick in the continuum over the 
entire time sequence, thereby ensuring at all times the adequacy 
of the diffusion approximation employed at this inner boundary. 
(Note that we refer to the photosphere as the location where the 
inward-integrated continuum optical depth equals two thirds.) To 
isolate the effects associated with time dependence, we neglect any 
chemical stratification (homogeneous ejecta) and assume a uniform 
power-law density distribution {d\np/d\nR = — n), with the ex- 
ponent n set to 10. 

We illustrate time-dependence effects with four different mod- 
els. We set up a "b aseline" model (A) that mimicks the proper- 
ties of SN 1999em jPessart & Hillieij|2006ah . In practice, we ini- 
tially map onto the CMFGEN grid an envelope with the follow- 
ing properties: Rq = 2.524 x 10^* cm (J?o is the base radius), 
vo = 3470 km s"^(uo is the base velocity), i?i^po = 6.781 x 10^^ g 
(po is the base density), p(R) = po (Ro/R)"' with the density ex- 
ponent n = W (R is the radius), and Lo — 3.6 x lO^^L0(Lo is 
the base luminosity). The maximum radius is 16 times the value of 
the base radius Ro (this choice ensures that at the outer CMFGEN 
boundary, the radiation field is more decoupled from matter at all 
frequencies and free streaming, i.e., the Eddington factor is close 
to unity). These initial conditions correspond to a base Rosseland 
optical depth of 3640, a total mass of 6.17Mq, and, assuming ho- 



mology, a time after explosion of 8.4 days. We adopt a fixed compo- 
sition (i.e,. chemical homogeneity), with values consistent with the 
red-sugergiant pro genitor of SN 1999em (Dessart & Hillier 2006^; 
ISmarttetal .1120021) . Adopted abundances, given by number, are: 
H/He = 5, C/He = 0.00017, N/He = 0.0068, 0/He = 0.0001, and 
all metal species are at the solar value. We also set up a slow (B; 
SN 1999br-like) and a fast (C) model, differing in the velocity of 
its inner mass shell, i.e., 2000km s^^(4400km s~^) rather than 
3471km s^^ The inner radius is also adjusted so that we start in 
these three cases at about one week after explosion. Finally, we 
set up an "early-time" model (D) that starts at ~2days after ex- 
plosion. In this model, we employ a slightly different composition 
for CNO elements, with C/He = 0.0001, N/He = 0.0013, 0/He = 
0.0016 (Dessart et al. 2007). We give a synopsis of all these model 
characteristics in Table[T] Apart from the "early" model that covers 
only the first week, we evolve these models from one week up to 
about 6 weeks after shock breakout. 

In our current approach, each ion-level population is evolved 
in a time-dependent non-LTE fashion, and thus, for simplicity, 
we adopt a unique model atom for all steps in a given model 
time sequence. For the baseline, slow, and fast models, we treat 
Hi (30,20), He I (51,40), Hell (5,5), ClI (59,32), Cm (20,12), 
CiV (14,9), N I (104,44), Nil (41,23), NlII (8,8), O I (75,23), O II 
(111,30), OlII (46,26), Nal (71,22), Mgll (65,22), Sill (59,31), 
Sim (53,29), Can (77,21), Till (152,37), Fell (309,116), Felll 
(477,61), FelV (282,50), and Nill (93,19), where the parenthe- 
sis contains the number of full- and super-levels (see Hillier & 
Miller 1998 for details). Compared to the model atom used in 
iDessart & Hillieil ( l2006al) . we omit Cl, neon, aluminium, sulfur, 
chromium, manganese, and cobalt. Given the much higher ion- 
ization state at earlier epochs, the "early" model neglects most 
low-ionization species (and metals other than iron) and adds He II 
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Figure 3. Top row: Compaiison of the radial variation of the ionization state of all species included in model A, for the time-dependent (left) and the steady- 
state case (right). Note how the effects on hydrogen and helium discussed in i|3.1l and illustrated in Fig. |2] apply in fact to all species, with a systematic 
over-ionization at large distances. Bottom row: Same as for the top-row panels, but this time shown in the log to reveal better the magnitude of the predicted 
difference in ejecta ionization between time-dependent and steady-state approaches. [See the electronic edition of the Journal for a color version of this figure, 
and the text for discussion.] 



(30, 13), C IV (64,59), N iv (60,34), N v (67,45), O iv (72,53), O v 
(152,75), Fe V (191,47), und Fe VI (433,44). 

From a converged initial steady-state model, we then step in 
time by Ai = 3 days (0.7 day for the early model D), i.e., all mass 
shells m in the initial model are evolved from their position R{m) 
to a new radius R'{m), with R'{m) — R{m) + vo{m)At, a new 
densityp' = p{R/ R')"^ , and the same velocity. In practice, for each 
new model, CMFGEN adapts its 80-point spatial grid to cover each 
optical-depth decade with at least five points. For model D, the grid 
contains only 50 radial points. 

In ij4] we present a sample of comparisons to observations to 
illustrate the performance of these time-dependent models (note 
that our models were not tailored to match any specific observa- 



tion). For the "early" model D, observations are only starting now 
to be available through, e.g., the SWIFT satellite (Brown et al. 
2007), and the corresponding simulations for model D, thus, re- 
main largely unconstrained (but see Quimby et al. 2007; Dessart 
et al. 2007). To illustrate the bolometric luminosity and color evo- 
lution of our four models, we present in Fig. [T] their light curves 
covering from the UV to the near-IR, with properties quite typi- 
cal of Type II-P SNe (Leonard et al. 2002; Pastorello et al. 2006; 
Dessart et al. 2007). 
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Figure 4. Radial variation of the electron density (solid line; in the 
logarithm, and in units of cm^'^) and the temperature (broken line; in units 
of 1000 K), for the time-dependent baseline model (black), and its steady- 
state equivalent (red). Both models have otherwise identical parameters and 
correspond to a time of 48.7 days after explosion. [See the electronic edition 
of the Journal for a color version of this figure, and the text for discussion.] 



3 RESULTS 

3.1 Time-dependence effects on ejecta properties 

We start our discussion of time-dependence effects on Type II SNe 
by focusing on the properties of the ejecta. In Fig. [2l we show 
the ionization fraction for hydrogen (left) and helium (right) for 
model A at the end of the time sequence, 48.7 days after the explo- 
sion. Hydrogen recombination occurs just above the photosphere 
but less efficiently in the time-dependent model than in the steady- 
state model, neutral hydrogen being 3 orders of magnitude more 
abundant in the later model. In the time-dependent model, two ad- 
ditional features of interest are the presence of fully ionized hydro- 
gen at large distances above the photosphere (the hydrogen ioniza- 
tion state is frozen) and the delayed recombination in the optically 
thick layers, at and below the photosphere (the ionization front is 
located further out). All these features apply to the helium ioniza- 
tion fraction as well, which we show in the right panel of Fig. [2] 
The Hell population is 15 orders of magnitude lower than that of 
Hel at intermediate heights, but it is comparable or stronger at large 
distances. Effects of time dependence can thus reach considerable 
levels, much beyond modulations that can be achieved through rea- 
sonable changes of (or uncertainties in) chemical abundance or ion- 
izing fluxes. 

Because metal line-blanketing plays an increasingly important 
role at and after the recombination epoch, time-dependence effects 
on metals are also of special interest. We show in Fig. [3] the ion- 
ization state of all species for model A on day 48.7. Notice how 
the frozen-in ionization is readily visible at large distances above 
the photosphere, with order of magnitude differences in ioniza- 
tion level between the time-dependent and steady-state approaches. 
Hence, the inclusion of time-dependence affects all species, irre- 
spective, for example, of the atomic properties of the given ion, and 
with a magnitude that is greatest at large distances above the pho- 
tosphere where the recombination timescale becomes increasingly 
large compared to the expansion timescale. 

Finally, we show in Fig. |4] the radial variation of the temper- 
ature and the electron density for model A at 48.7 days after ex- 
plosion. The temperature in the time-dependent model is higher 
in the optically-thick regions, due to the reduction of density- 



sensitive line emission processes that dominate radiative cooling. 
In the optically-thin regions, above the photosphere, the tempera- 
ture is, however, quite similar between the two approaches. This 
result stems from the poor coupling between gas and radiation in 
these scattering-dominated, tenuous, ejecta. Departures from LTE 
are pronounced and the temperature is a secondary driver of the 
ejecta ionization, even in steady-state configurations. 

Models B-C-D offer valuable information on the sensitivity 
of the above effects on the ejecta velocity (models B and C) and 
the epoch (model D). We find that varying the velocity of the 
ejecta within the range of values for slow/underluminous (model 
B; SN 1999br) and more standard luminosity Type II SN events 
(SN 1999em) has no sizeable effect on the main features of time 
dependence highlighted above. Recombination is inhibited, and a 
significant over-ionization persists far above the photosphere. 

In Fig.|5] we present the mean ionic charge for the six species 
(H, He, C, N, O, and Fe) included in model D (early epoch), and at 
6.7 days after explosion. While hydrogen remains fully ionized, in 
both time-dependent and steady-state approaches, all other species 
show a higher ionization at large distances above the photosphere; 
no difference is noticeable in the vicinity and below the photo- 
sphere, contrary to model A at 48.7 days. This supports the idea that 
time-dependence effects operate even at early times. They do so 
at large distances because there, the recombination and expansion 
timescales are already sufficiently close, but they can also act at the 
photosphere because of optical depth effects and meta-stable lev- 
els (see i]2.3t . It also shows that over-ionization at large distances 
does not stem by essence from the recombination energy (since 
most of that energy is stored in hydrogen and hydrogen remains 
fully ionized in this whole sequence) but results exclusively from 
the low-density fast-expansion of the SN ejecta. A corollary is that 
this progressive "over-ionization" does not require an extra source 
of ionization and excitation from, for example, ^®Ni (Mitchell et al. 
2001; Baron et al. 2003), and operates instead, in a robust fashion, 
through the reduction of the effective recombination rates. Finally, 
this supports the idea that similar effects should apply in SN ejecta 
in general, since they generically possess similar low density and 
large velocities. 

We conclude this section by presenting the evolution of the 
electron-scattering optical depth at the base radius Rq of each 
model in Fig. |6] (see Table [T] for model parameters). Rather than 
having a dependence (dashed lines) for an expanding ejecta 
with a fixed mass absorption coefficient (as in model D where the 
dashed and solid lines overlap; fixed composition and ionization), 
the curves for models A, B, and C show a change of concavity, 
with a slow decrease at first and a turn-over to a faster decrease 
after 3-6 weeks. Recombination, which is delayed at first and per- 
mits the persistence of a high free-electron density in the ejecta, 
does in the longer term occur and triggers a fast drop of the base 
electron-scattering optical depth. It remains, however, higher by up 
to a factor of about ten, than the value obtained in the correspond- 
ing steady-state model (for the same adopted lower bound radius 
in the optical-depth integral) at the ultimate time in the sequence, 
shown as filled colored circles. Model D, whose ejecta contains 
species at least once ionized throughout the sequence, has a base 
optical depth of 51.9 in the time-dependent approach and 51.7 in 
the steady-state approach, the modest difference stemming from 
the slight over-ionization of sub-dominant species like helium and 
CNO elements. (See the last columns of Table[T]for a log of these 
values.) 
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Figure 5. Comparison of the radial variation of the ionization state of all species included in model D (the model that covers the first week of evolution of 
a Type II SN ejecta), for the time-dependent (left) and the steady-state case (right). Note how the behaviour illustrated in Fig.|2]for hydrogen and helium is 
visible in fact in all species (with the exception of hydrogen which remains fully ionized in both time-dependent and steady-state approaches at such an early 
time; black line). Hence, even during the first week after explosion, when the ejecta density are relatively higher, time-dependence effects leave an obvious 
imprint on the ejecta ionization state (He, C, N, O, Fe here), even though the medium, made-up primarily of fully-ionized hydrogen, retains roughly the same 
electron density. [See the electronic edition of the Journal for a color version of this figure, and the text for discussion.] 
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Figure 6. Time evolution of the electron-scattering optical depth To.s. 
for each time-dependent model A (black; baseline model), B (blue; slow 
model), C (green; fast model), and D (red; early epoch model), using sohd 
lines and cross symbols (to mark the days computed). We oveiplot the elec- 
tron scattering optical depth for the corresponding steady-state model at the 
final time in each sequence (filled circle). For each model, we also draw 
the 1/i^ optical-depth evolution for an ejecta with fixed electron-scattering 
opacity (fixed composition) and ionization (dashed line; the blue and black 
dashed curves overlap, and so do the red-dashed red-solid curves). [See the 
electronic edition of the Journal for a color version of this figure, and the 
text for discussion.] 



3.2 Time-dependence effects on the emergent light, from UV 
to near-IR, for the baseline model A 

The various investigations and analyses discussed above indicate 
that time-dependence effects are robust, apply to all species, even 
prior to the recombination epoch of hydrogen. Hence, for example, 



during the first week after explosion, such effects inhibit the re- 
combination from Fe''"^ to Fe^"^ and Fe'^"'", or that of He^ to He""*". 
Different species block light at specific wavelengths, with often a 
tendency to bunch-up in specific spectral regions, but with no a- 
priori correlation, and therefore we anticipate that such modula- 
tions in the ionization state of the ejecta might lead to a significant 
impact on the properties of the emergent light. Both the optical, 
dominated by Hi and Hel lines, and the UV, dominated by metal 
line-blanketing due to, for example, iron, should be affected. 

In Fig. [7] we present a montage of synthetic spectra showing 
the time evolution for the time-dependent (black) and steady-state 
(red) CMFGEN model A in the UV/optical (left panel), the 9000 A- 
-2pm region (middle panel), and the 2-5pim region (right panel) at 
three consecutive epochs (8.41 d, the start of the sequence, 30.6 d, 
and the last time in the sequence at 48.6 d). The SED has a peak that 
shifts from the UV to the optical over the 6 weeks covered, starting 
first as a nearly featureless continuum and then revealing more and 
more the effects of metal line blanketing in the UV and optical (see, 
e.g.. Brown et al. 2007). At a qualitative level, this model sequence, 
which accounts for time dependence, is in agreement with what 
is observed in nature, for example for the Type Il-P SN 2006bp 
(Dessart et al. 2007), and what could also be modeled with steady- 
state CMFGEN models. Hence, the time-dependence incorporated 
here has no major influence on the general morphology of the SED 
and its evolution with time. A corollary is that time-dependence 
effects on the continuum SED are weak. The effects on lines are, 
however, pronounced as we notice that all hydrogen lines, from the 
Balmer to the Paschen and to the Brackett series, are considerably 
weaker in the steady-state model for the last date, compared to its 
time-dependent equivalent, suggesting that effects do not intervene 
selectively, but rather that all hydrogen lines are affected. 

To delineate more precisely the influence on lines, we display 
a montage of optical synthetic spectra for model A, correspond- 
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Figure 7. Montage of (scaled) synthetic spectra showing the time evolution for the time-dependent (black) and steady-state (red) CMFGEN model A in 
the UV/optical (left panel), the 1-2/im region (middle panel), and the 2-5^m region (right panel). A vertical shift is applied to each spectrum, smaller with 
increasing time. We show the evolution for three dates, at 8.41 d (top curve; the time-dependent and the steady-state both start their evolution at that time and, 
thus, only one curve is plotted), slightly past halfway through the time sequence at 30.6 d (middle curve), and for the last time computed at 48.6 d (bottom). 
[See the electronic edition of the Journal for a color version of this figure, and the text for discussion.] 



ing to the last time in the sequence, in Fig. [S] In the top panel, 
we plot the result for the time-dependent (black) and the equiva- 
lent steady-state (red) model, while in the lower panel, we present 
rectified spectra obtained by post-processing the CMFGEN results 
and accounting for all continuum processes but only bound-bound 
transitions of selected ions. The contrast for Hydrogen Balmer lines 
is stark, and most extreme for Ha, which is now as strong as ob- 
served. Helium, which does not show a single line in the steady- 
state model, shows a non-trivial absorption feature at ~1.03/^m 
and a very extended weak flat-topped emission profile (we will 
return to this line below). Continuing upwards, nitrogen and oxy- 
gen reveal only weak lines, and effects are marginal. Steady-state 
models show stronger absorption and weaker emission in Nl lines. 
The same holds for Mgll and Sill. In the time-dependent case, 
Nal5895A strengthens and broadens. This occurs because Na"*" is 
the dominant ionization state in both models, but the increased elec- 
tron density in the time-dependent model leads to an increase in the 
neutral Na fraction. Conversely, Call 8500A weakens, primarily be- 
cause the relative fractions of Ca^ and Ca"''^ change significantly. 
Fell and Till lines are stronger in the steady state model — stronger 
in absorption in the UV and stronger in emission in the optical — 
resulting from the lower ionization of the ejecta. When combined, 
all these effects give an SED that is most conspicuously different 
through the changes in the strength of Ha and Call 8500A, whose 
ratio is reversed between the two approaches. 

An alternative approach of strengthening the weak Ha pro- 
file, through a flattening of the density distributiorQ, leads instead, 
primarily, to a considerable strengthening of the Call 8500A multi- 
plet, in complete contradiction with observations jPessart & Hillierl 
l2005ah . Time dependence offers here a means to produce a strong 
Ha line, while delayed recombination of calcium maintains the 
moderate Call 8500A strength. The moderate impact on metals in 

^ Note that the recombination efficiency scales as the product of the elec- 
tron density with the ion density rii. Flattening the density distribution 
enhances the material density p above the photosphere, but if the ejecta is 
recombined, the low values of rie and rii will inhibit the recombination 
needed to power, for example, a strong Ha line (see also S|4). 



the vicinity of the photosphere leaves metal line-blanketing only 
weakly affected, and it is only the strongest lines (like Fell 5 169 A) 
that show both a strengthening and a broadening in the less-ionized, 
steady-state, model. 

In Fig.[9l we plot the radial distribution of the parameter C, of 

iHillieil JI987D for the Ha line, a quantity that relates to the total 
emission in the line through the integral J^"''"' C,{R')d\og R' (or, 

equivalently, jj"" ^(V^')dlog V' for homologously expanding 
ejecta). The curves for the time-dependent (solid line, in black apart 
from the final model in red) model A indicate a sustained and com- 
parable emission distribution throughout the sequence. The peak 
emission shifts from 8000km s^^at one week to 4000km s^^at 6 
weeks after explosion, following the expansion of the ejecta and the 
decrease in density. But the emission region considerably broad- 
ens with time, so that, while its inner edge recedes, its outer edge 
remains fairly fixed (in a Lagrangean sense). The corresponding 
emission region for the steady-state model (broken red line) is lo- 
cated deeper (the photosphere as well; see i]3.1l > and is considerably 
narrower. This property translates directly in the observed Ha line 
properties, with a smaller line flux (weak line) confined to smaller 
velocity (narrow line) in the steady- state model. 

To f urther illustra te the increase of emission-line volumes, and 
following iDessart & H illier (2005a, b), we show in Fig.[TOl for H7 
(left-column panels), H/3 (middle-column panels), and Ha (right- 
column panels), grayscale images in the {v,p) plane of the flux-like 
quantity p ■ 7(p), where II = [(A/Aq) — l]c is the classical Doppler 
velocity, p is the impact parameter in units of the photosphere ra- 
dius -Rphot, and I{jp) is the emergent specific intensity along p. 
The sum over p of the quantity p ■ I{p) aiv corresponds to the to- 
tal line flux at v, shown at the top of each panel (solid line). Note 
that for each line, we select a single Hi transition to avoid the cor- 
rupting effect of line overlap as we near the Balmer jump (all other 
bound-bound transitions of Hi and of other ions a re ignored). A de- 
tailed in terpretation of these diagrams is given in lPessart & Hilliea 
( l2005j|bl) and will not be repeated here. Note the stark contrast 
between the considerably more extended regions of emission and 
absorption (in velocity or equivalently in space dimension) in the 
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Figure 8. Top: Comparison of synthetic spectra obtained for model A and 
the last time in the sequence at 48.7 days, showing the results for the time- 
dependent model (black curve) and the equivalent steady-state model on the 
same day (red curve). Bottom: Same as top, but for the synthetic rectified 
spectra (shifted vertically for visibility, and magnified by a factor of five 
beyond 5500A for all species labelled with a superscript a) and allowing 
exclusively for bound-bound transitions of the ion species indicated on the 
right. We only show those species that leave an imprint on the spectrum. 
This figure illustrates the modest time-dependence effects over the first 6 
weeks after explosion on the continuum. However, the effects on lines are 
significant, most noticeable on the strengthened hydrogen Balmer lines and 
NalD, and the weakened Call 8500A multiplet lines. The enhanced recom- 
bination of the ejecta in the steady-state model results in a lower ionization 
state and stronger metal lines, e.g., for Fell and Till. [See the electronic 
edition of the Journal for a color version of this figure, and the text for dis- 
cussion.] 



time-dependent models (top row) compared to the very-confined 
absorption and nearly nonexistant emission for the corresponding 
lines in the steady-state counterparts (bottom row). The same ejecta 
is covered by both models (same mass shells with same velocities), 
yet the Balmer lines show considerably broader profile in the time- 
dependent case, mimicking an ejecta with a faster expansion rate at 
that given time. 

The time-dependent version of CMFGEN, thus, resolves the 
discrepancies with hydrogen Balmer lines obtained with its steady- 
state version. While the difficulty with Balmer lines has been docu- 
mented in the past (see, e.g., UC05), we find that time-dependence 
effects conspire to modify the strength of most lines, not just those 
of hydrogen. For example, Nal D undergoes a considerable broad- 
ening and strengthening with the treatment of time dependence, 
while in the past, fitting it well r equired enhancing the so dium 
abundance by a factor of about four jPessart & Hiniej2006ah . Per- 
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Figure 9. Solid line: Time evolution , in velocity space and for model A, 
of the 5 parameter of iHillieJ i 19871) for Ha, which describes the spatial 
distribution of the corresponding line emission (the area under the curve 
scales with the emergent line flux). Notice how this distribution peaks at 
smaller radii for later times, but also the persistence of emission at large 
radii at all times (the last time shown in red). Brolien line: Same as for 
the solid line but for the steady-state CMFGEN model, and only for the 
last time (to be compared to the solid red line). Note the deeper location 
of the emission, resulting in a very narrow Hq profile, in contradiction to 
observations at the conesponding mid-photospheric-phase epoch. [See the 
electronic edition of the Journal for a color version of this figure, and the 
text for discussion.] 

haps the most extreme case in our investigation is Hel 10830A, 
which persists over the entire 6 weeks covered by the baseline 
model A, and for which, at the last time in the sequence, we predict 
a non-trivial absorption strongly blueshifted from line center (at 
1.038/im, equivalent to -13000 km s^^), and a weak and broad flat- 
topped emissior0. Such a line profile morphology is synonymous 
with absorption/emission far above and detached from the photo- 
sphere, and usually associated with chemical stratification of the 
corresponding element or density kinks; here it stems solely from 
an ionization stratification caused by time-dependence effects. The 
Hel 10830A line is entirely absent in the steady-state counterpart 
model three weeks after the explosion because, in the relevant re- 
gion, helium atoms are once ionized with a probability of only one 
in lO^"-*. 

In Fig. [TT] we reproduce Fig. [TO] and illustrate the absorp- 
tion/emission sites for Hel 10830A at 48.7 days in the time- 
dependent version of IVIodel A. Absorption stems from the re- 
gions in the direction of the stellar disk, far from the photosphere 
which has a velocity of ~4500km s^^on that day. For impact 
rays not intersecting the photodisk, only emission occurs, again at 
large distances above the photosphere, and over a small volume. 
The emission flux is small, appears as a flat-topped profile, and is 
thus barely noticeable above the continuum. Of this Hel 10830A 
line in the synthetic spectrum, one can only observe the strongly 
blueshifted absorption component, the location of the flat-topped 
emission overlapping with the much stronger P7 line at this late 
epoch. This is illustrated differently in Fig. [121 where we show 

^ Chugai et al. (2007) propose that the cool dense shell that forms at the 
interface between the SN ejecta and the pre-SN wind is at the origin of 
a similar, but stronge r, feature in Hel 10830A, that fits more sui tably the 
observed line profile j Hamuv et al.ll200ltlSpvromiho et alj|l99lh . We ac- 
knowledge this possibility, and focus here oiily on the time-dependence ef- 
fects germane to the ejecta alone. 
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Figure 11. Same as in Fig. [To] but for the Hel 1O830A line (to enhance 
the rendering, we subtract p ■ I(p) from the continuum-only equivalent, 
and add a few contours). The line profile is flat-topped, shows a "de- 
tached" blueshifted absorption, forms further out than Ha in the same time- 
dependent model, and is completely absent in the equivalent steady-state 
model. [See the electronic edition of the Journal for a color version of this 
figure, and the text for discussion.] 

a montage of synthetic spectra covering the 1-1.2/im region and 
the entire evolution computed for Model A. In the left panel, we 
show the full synthetic spectrum, i.e., including all species. In the 
right panel, we use the converged CMFGEN model at each time in 
the sequence and solve the formal solution of the radiative trans- 
fer problem by including bound-bound transitions of helium only. 
Note how at early times (dark curves) the profiles in the two pan- 
els match, indicating that helium is the main feature contributing 
in emission and absorption at early times. After three weeks, the 
emission is dominated by P7 and the Hel absorption gives rise to 
a dip, detached from the P7 absorption, and shifting further to the 
blue (i.e., in an opposite sense to the photosphere which recedes to 
smaller velocities, as well as to all other line-profile absorptions). 

3.3 Sensitivity to ejecta properties and epochs 

Besides the baseline model A discussed above, we have computed 
three other time sequences, as described in detail in ^2.5\ Model B 
and C investigate the effect of the expansion velocity of the ejecta, 
model B being also an underluminous event by design. Model D 
further investigates the possible (but not expected) time-dependent 
effects during the first week after explosion. 

We show in Fig. [T3] a comparison between optical synthetic 
spectra computed for the time-dependent (black) and the steady- 
state (red) Model B (left) and Model C (right), both at the last time 



in the sequence. Despite their different expansion velocities, mod- 
els B and C show very similar time-dependence effects, with the 
most noticeable features being the same as for model A over the 
same time range: Hydrogen Balmer lines are strengthened, and so 
are those of the Paschen and the Brackett series (not shown here 
but see Fig. [7] for a guide), as well as Nal D. Here again, and 
by contrast, the Call multiplet at ~8500A is weakened. Impor- 
tantly, given that these models correspond to the conditions from 
underluminous to standard, and from slow to fast ejecta, germane 
to SNe 1999br/2005cs and 1999em/2006bp, time-dependent CMF- 
GEN models systematically predict strong spectroscopic effects. In 
an absolute sense, these modulations still correspond to ejecta that 
are fast and tenuous, where the expansion time scale is comparable 
to the magnitude of important rates (like that of recombination), 
affecting strongly the ionization equilibrium at and above the pho- 
tosphere. 

Turning to Model D in Fig. [14] we see that the differences, at 
6.7 days, between the time-dependent and steady-state model SEDs 
are subtle (bottom panel). The continuum SED is essentially iden- 
tical and noticeable differences are only visible in the Hel lines, 
particularly at 5875A and at 10830A. Both of these are a factor of 
about 2-3 stronger in the time-dependent model. In the past, repro- 
ducing the strength of optical Hel lines with a typical BRG/RSG 
surface helium abundance has been a challenge (Eastman & Kirsh- 
ner 1989), although this difficulty disappears with a full non-LTE 
treatment (Dess art & Hilliej2005ah . We now find that time depen- 
dence introduces a further tuning, making optical Hel lines persist 
over a longer time span than is predicted by (non-LTE) steady-state 
models. 



4 COMPARISON WITH A SAMPLE OF OBSERVATIONS 

The motivation for the present work was to find an explanation, 
and a cure, for our inability with steady-state CMFGEN models 
to reproduce the strong Ha line at the recombination epoch. De- 
spite all attempts, we were unable to reproduce the Ha line strength 
as soo n as hydrogen recombina tion in the Type II SN ejecta takes 
place jPessart & Hillieill2b06ah . We reported this problem for SN 
1999em, but it is in fact a generic feature of all steady-state CM- 
FGEN models at a few weeks after explosion. Interestingly, obser- 
vations at the corresponding epoch reveal that Ha is instead the 
strongest line in the optical spectrum of Type II SNe, across the 
whole range of objects, from the peculiar SN 1987 A, to the un- 
derluminous slowly-expanding SN 1999br, to the more standard 
SN 1999em (Type II-P). These SNe are understood to stem from 
core-collapse explosions of supergiant stars, and are characterized 
by a wide range of explosion energies (velocity distribution versus 
Lagrangean mass), nickel yields, and bolometric luminosity. The 
mechanism that can solve the present Ha discrepancy must operate 
at the hydrogen recombination epoch and be effective in all Type II 
SN ejecta. 

Specifically, it is unlikely that, in general, non-thermal ion- 
ization/excitatio n by ^^Ni is the missing mechanism in our for- 
mer approach ( [Dessart & Hillieil l2006al) since SNe events like 
1999br or 2005cs, which both boast a strong Ha line at the hy- 
drogen recombination epoch, are also characterized by vanishingly 
small ^"^Ni yields (Pastorglloet al. 2004; Tsvetkov et al. 20o|; 
iTakats & Vinkdll20oA|Pastorello et al. 2006). Furthermore, the Ha 
discrepancy appears at times when the photosphere still moves at 
4000-5000 km s~^, vel ocities out to which^^Ni is not expected 
to be mixed routinely dFrvxell et all 1 199 ll : iKifonidis et al.l l2000l 
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Figure 10. Grayscale images for H7 (left column), H/3 (middle column), and Ha (right column), of the quantity p ■ I{p) as a function of impact parameter p 
and classical Doppler velocity ti = (A/Ao — l)c, where p is the impact parameter (in units of the photospheric radius Rphot). ™d I{p) is the specific intensity 
along p (at v). The panels in the top (bottom) row correspond to the time-dependent (steady-state) results obtained for the model A ejecta at 48.7 days after 
shock breakout. The top half of each panel gives the line profile flux, which, for each v, corresponds to sum of p ■ I{p) over all ps, giving a vivid illustration 
of the (p, v) sites contributing to the observed line profile. To emphasize the differences between lines and model assumptions, axis extrema are kept identical 
for all frames. 



l2003h . The Ha problem in SN 1987A, for example, is already visi- 
ble at just a few days after explosion, while Cobalt lines, testifying 
for the presence of radioactive material in the photosphere layers, 
were first seen at about 45 days after explosion (see, e.g.. Pinto 
& Woosley 1988), significantly later (but much earlier than antici- 
pated by the theory at the time SN 1987 A went off). 

We start the comparison with optical and near-IR observa- 
tions of SN 1999em on the 18th- 19th of November (Leonard et 
al. 2002, Hamuy et al. 200 1). These observations were modeled 
bv lDessart & HillieJ ( 120063) and we reproduce their Fig. 7 in the 



right panel of Fig. [15] In the left panel, we fit the observations with 
the synthetic spectrum computed for model A on day 35, and in- 
clude the near-lR range to emphasize the concomitant strength- 
ening of the Paschen lines in the time-dependent model. Time- 
dependence effects are thus not limited to Ha, nor to hydrogen 
B aimer lines exclusively, but operate as well on the Paschen se- 
ries, and, even more generally, affect lines of all species. Note also 
that the predicted Hel line at 10830A is too weak to fit the ab- 
sorption featureatj_^3/jm, which, despite our neglect of Cl contri- 
butions jPessart & Hillier 2006a) , may require the additional con- 



14 Luc Dessart and D. John Hillier 




10000 10500 



11000 
A (A) 



11500 12000 



— I E 1 ] 3 1 1 1 1 1— 

Days oflSf Explosion 




10000 10500 



1 1000 
A (A) 



11500 12000 



Figure 12. Left: Time sequence of the rectified syntlietic flux over tlie 10000-12000A range including all the bound-bound transitions. Main features are 
due to P7 and Hel 10830A. Right: Same as left, but this time, we only include Hel bound-bound transitions. Notice the flat-top and blueshifted absorption of 
Hel IO830A at late times. [See the electronic edition of the Journal for a color version of this figure, and the text for discussion.] 
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Figure 13. Optical synthetic spectra for the time-dependent (black) and steady-state (red) model B (left; day 42.8) and model C (right; day 35.5) at the last 
time in the sequence. [See the electronic edition of the Journal for a color version of this figure, and the text for discussion.] 



tributio n from the Co ld Dense Shell at the ejecta/pre-SN-wind in- 
terface jChugai et alj.2007.) . Note that the time-dependent model, 
unlike the steady-state model shown on the right, was not tuned 
to fit the observations. Hence, the adopted luminosity evolution 
may not be fully adequate for this event. An additional issue is 
that model A started at one week after explosion, while we demon- 
strated clearly in i]3.2l the strong time-dependence effects occuring 
for Hel even at such early times, effects that will carry over to sub- 
sequent times and that may therefore be underestimated in these 
model sequences. 

In Fig. [16] we compare synthetic spectra for IVIodel A at 48.7 
days after explosion, with allowance for time dependence (left 
panel) or assuming steady-state (right panel), with the observa- 
tions of SN 1999em taken on the 15th of December 1999 (Leonard 
et al. 2002). The SED has changed significantly from that seen 
on the 19th of November 1999. There is more flux depletion in 
the UV, more metal line-blanketing in the optical, and Ha, Nal D, 



and Call 8500A have strengthened. The steady-state model (right 
panel) reproduces well the overall shape of the SED, but what was 
a slight discrepancy in fitting the Ha line at earlier epochs now be- 
comes a major discrepancy. By contrast, the time-dependent model 
reveals a strong and broad Ha line that closely matches the ob- 
served profile. The same agreement is reached for NalD. A num- 
ber of metal species are missing (scandium, barium) and the CNO 
abundances would also need adjustment (enhancement of carbon 
and oxygen abundances, and depletion of nitrogen; Dessart et al. 
2007) but overall, the time-dependent model allows a much more 
satisfactory fit to observations and in particular Ha, which is the 
defining feature of Type II SN spectra. 

Note that the alternatives, an enhancement in turbulent ve- 
locity or flattening of the density distribution, do not improve the 
agreement with observations. The former has a weak or negligi- 
ble effect at late times, while the latter leads to a large enhance- 
ment of the Call 8500A line strength, in stark contradiction with 
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Figure 15. Left: Illustration of the match between observations (black) of SN 1999em on the 19th of November 1999 in the optical range (top; Leonard et al. 
2002) and on the 18th of November 1999 in the near-IR range (bottom; Hamuy et al. 2001), with the reddened synthetic spectrum (red; blue: continuum-only) 
for the time-dependent model A at 35 days after explosion. We use the extinction law of Cardelli et al. (1988), with E(B — V) = 0.1, and normalize the 
synthetic flux to the observed flux at 6IOOA, with a fuilher 15% flux reduction in the near-IR (note the one day time difference between the two spectra). 
Right: Same as left, but this time using the equivalent steady-state model. [See the electronic edition of the Journal for a color version of this figure, and the 
text for discussion; see also §3.2.7 and Fig. 7 in Dessart & Hillier (2006a).] 



observations. The key issue here is that one should not view the 
recombination rate as scaling with the density squared, but rather 
as the ion density times the free-electron density. Once the mate- 
rial has fully recombined, the absence of free-electrons and ions 
completely shrinks the recombination rate, and Ha becomes only 
weakly sensitive to modulations in the outer mass density distribu- 
tion. 



We finally show in Fig.[T7]a comparison of Model B at 42.8 
days after explosion, including time dependence (left panel) or 
assuming steady state (right panel), with the observations of SN 
1999br on the 19th of May 1999 (Pastorello et al. 2004). The time- 
dependence effects are similar to those described for SN 1999em on 
the 15th of December 1999, but they bear more strongly on the role 
of time dependence since, for the underluminous low-nickel-yield 
SN 1999br, non- thermal ionization/excitation is very unlikely to af- 
fect in any way Ha emission. Here again, time dependence effects 
lead to strong Balmer (and Paschen) lines, as observed, without 
recourse to extra ionizing/exciting sources or changes in the den- 
sity profile. Furthermore, this shows that time-dependent effects are 
present in a wide variety of SN ejecta, even those characterized by 
a slow expansion rate. 



5 IMPLICATIONS 

5.1 Abundance determinations 

The persistence of Hel 10830A at 6 weeks after explosion in Model 
A suggests that the so-called detached absorption/emission lines 
are not exclusively associated with a stratification of abundances, 
with specific ions present at certain heights or mass shells in the 
ejecta. Here, the detached absorption/emission from Hel 10830A 
stems purely from a time-dependence effect, and thus does not in 
principle require chemical stratification. Time-dependence effects 
can alter our interpretation of chemical stratification based on line 
profile measurements, e.g., velocity at maximum absorption. This 
is supported further by the time-dependence effects we see on the 
profile shape of strong lines, such as, for example, the hydrogen 
Balmer lines or Nal D. In the time-dependent model A at 48.7 days, 
we measure from the synthetic spectrum that included only the 
bound-bound transition of the relevant species a velocity at max- 
imum absorption of -7500km s~^ for Ha, -5800km s~^ for H/9, 
-4900km s"^ for H7, -4000km s"^ for NalD (at 5890A), while 
in the same order, measurements performed on the synthetic spec- 
trum for the steady-state model give -3900km s^^, -3500km s^^, 
-3400km s~^, and -3200km s~^. Differences between these val- 
ues are greater than those in photospheric velocity between the two 
models, which, respectively, are 4130km s~^ and 3510km s"\ 
Time-dependence effects operate in a non-linear fashion, affect- 
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Figure 14. Top: Comparison of synthetic spectra obtained for model D at 
6.8 days for the time-dependent (black) and steady-state (red) approaches. 
Bottom: Same as top, but for the synthetic rectified spectra accounting only 
for Hel bound-bound transitions. No time-dependent effect in the emergent 
light is seen for the other species treated in the calculation. (See text for 
discussion.) 



ing lines and the continuum with a different magnitude, something 
to be expected since the continuum is mostly sensitive to the ion- 
ization state (a global quantity) while line profiles are affected by 
many parameters, such as the relevant ion level populations, the 
electron density, the radiation field etc. over a broad region in the 
SN ejecta. 

Tomographic techniques, which aim at constructing the distri- 
bution of elements with mass in SN ejecta through the assumption 
of homologous expansion and kinematic measurements on (or fits 
to) line profiles (Stehle et al. 2005), may be compromised through 
the neglect, in the solution to the radiative transfer problem, of 
time-dependent terms in the statistical and radiative equilibrium 
equations. In the future, we will investigate such effects on infer- 
ences of chemical stratification in SN ejecta. 



5.2 Correction factors 

Type II SNe offer an attractive means to determine distances to 
galaxy hosts, and ultimately, constrain distances in the Universe. 
The so-called Expanding Photosphere Method (EPM; Kirshner & 
Kwan 1974) has been used and improved upon to yield distances 
to a few nearby Type II SNe (Eastman & Kirshner 1989; Schmidt 
et al. 1994; Eastman et al. 1996), although with a moderate level 
of success for SN 1999em for which the EPM distance (Hamuy 
et al. 2001; Leonard et al. 2002a; Elmhamdi et al. 2003) under- 
estimates by 50% the Cepheid distance (Leonard et al. 2003) to 
the galaxy host. Using the updated correction factors of Dessart 
& Hillier (2005b) partially cures this discrepancy, while even more 
consistent results require the use of tailored models, with their com- 
bined correction factors and color temperatures (Dessart & Hillier 
2006; Dessart et al. 2007). 

Originally, the correction factors were introduced to account 
for the flux dilution that stems from the large electron-scattering 
optical depth contribution of the photosphere. However, as noted 
by Dessart & Hillier (2005b), they also correct for the (strong) de- 
viations from a blackbody distribution arising from the UV and 
optical line blanketing, particularly at later times. 

Here, we investigate whether time-dependent models yield 
correction factors and associated color temperatures that depart 
systematically from steady-state models. Time dependence affects 
flux dilution by enhancing the electron density (and thus the magni- 
tude of flux dilution) at all times, relative to a steady-state approach. 
However, it also affects the SED directly, both for the continuum 
(weakly) and for lines (strongly). We show in Fig. [TS] the collec- 
tion of correction factors versus color temperature for each set BV 
(left), BVI (middle), and VI (right), for all models computed in 
the present study (51 models, shown as dots). For comparison, we 
overplot correction factors and associated color temperatures com- 
puted with steady-state CMFGEN models (37 in total, shown as 
crosses) for the study of SN 1999em (Dessart & Hillier 2006a}, as 
well as that for SNe 2005cs and 2006bp (D essart et al.ll2007h . 

At temperatures above 10,000K we find no sizable systematic 
difference in correction factors computed with the time-dependent 
and steady-state approaches, which supports the use of the steady- 
state models of Type II SNe for distance determinations. At lower 
temperatures, particularly at temperatures less than 7000K, the new 
correction factors are somewhat smaller, and tend to lie at the 
lower edge of the distribution of correction factors calculated from 
steady-state models. They are also closer to those of Hamuy et al. 
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Figure 17. Same as Fig. 1151 but for the observations of the underluminous event SN 1999br (black line; Pastorello et al. 2004) and using the time-dependent 
model B at 42.8 days after explosion. Synthetic spectra are reddened using the extinction law of Cardelli et al. (1988) with E{B — V) = 0.03. [See the 
electronic edition of the Journal for a color version of this figure, and the text for discussion.] 



(2001). The large scatter from model to model suggests that a tai- 
lored analysis is nonetheless preferred over the use of analytical 
correction factors. The departures we observe at low temperatures 
are not surprising — this is a regime where the time-dependent 
terms have a very large effect on spectral features, and a regime 
which we have tried to avoid in our EPM distance determinations. 
The inclusion of the time-dependent terms will allow us to better 
model the long-term evolution of SN spectra, and hence to extend 
the time baseline over which the EPM can be used. Further, the in- 
clusion of these time-dependent terms is also important if we are 
to reduce systematic errors, and realize the fullest potential of the 
EPM technique. 



5.3 Polarization 

Through its effect on the ionization equilibrium of Type II SN 
ejecta, time dependence leads to an enhancement in the ejecta elec- 
tron density. The ejecta optical depth falls off more rapidly than 
but not as fast as it would for a steady-state configuration 
(Fig.[6l(. Deeper layers should thus be revealed with a delay and the 
nebular phase be entered later in nature compared to steady-state 
predictions. 

Moreover, the approximate linear scaling of the de- 

i ;ree of polarization wi th the electron-scattering optical depth 
Brown & McLeanll977h suggests that, for a given degree of asym- 
metry, higher polarization measures will result from theoretical 
models generated assuming time-dependence, rather than station- 
arity. In Type II-P SNe, polarization is usually inferred to be at the 
0.1% level during the plateau phase, and, although late time po- 
larization measures are rare, SN 2004dj indicated that polarization 
can increase dr amatically (here to ^--^1%) when the core is revealed 
( [Leonard et alj|2006,) . Understanding what controls the ionization 
at such late times, in particular the respective roles of time depen- 
dence and radioactive decay from ^^Co, is therefore critical for a 
proper interpretation of the degree of polarization from SN obser- 
vations and the link to the asymmetry of the explosion. 

6 CONCLUSION 

We have presented an extension to the non-LTE model atmosphere 
code CMFGEN to treat time-dependent terms in the statistical and 



radiative equilibrium equations. We use implicit first order differ- 
encing, and the resulting equations are solved by a partial lineariza- 
tion technique in a manne r virtually identical to that used to solve 
the steady state equations ( lHillier&MiTietlll998h . We confirm the 
findings of Utrobin & Chugai (2005), which were related to Ha 
and Ball 6I42A at a few days after explosion in the spectrum of SN 
I987A, that the time dependent terms are important for the analysis 
of Type II SN spectra. 

We find that the inclusion of time-dependence terms in the 
statistical equilibrium equations produces ejecta, in Type II SN, 
that are systematically over-ionized relative to steady-state mod- 
els. Spectroscopically, the associated changes in level populations 
and optical depths alter the line profiles of all species throughout 
the spectral range. Qualitatively, lines appear in general stronger 
and broader, as is seen for the Balmer, the Paschen, and the Brack- 
ett series of hydrogen or for Nal D. Exceptions (e.g. Call 8500A) 
arise when the altered ionization inhibits the necessary recombina- 
tion of a given ion (e.g., Ca^"*"). Because of optical depth effects, 
the inclusion of time-dependent terms in the statistical equilibrium 
equations can even influence the ionization state of the gas in re- 
gions where the classic recombination times-scale is much shorter 
than the flow time-scale. Surprisingly, we predict the presence of 
Hel 10830A 6 weeks after the explosion, naturally caused by the 
frozen, and full, ionization of the outer ejecta. The resulting de- 
tached blueshifted absorption and flat-topped emission profile is 
supported by observations, although the magnitude of the absorp- 
tion component may require the additional contribution from the 
cold dense shell that c an form at the interf ace between the SN ejecta 
and the pre-SN wind IChugai et alj2007h . 

The effects of time dependence observed stem from the energy 
gain that follows changes in ionization and excitation, and from the 
similarity in recombination and expansion time-scales (crucial at 
large distances above the photosphere). The recombination time- 
scale can be lengthened by optical-depth effects, and in some cases 
this may allow time-dependence effects to be important even at the 
photosphere. The present study indicates that neglecting time de- 
pendence may compromise the results of quantitative analyses of 
Type II SN spectra, affecting the inferred ejecta ionization (and its 
origin) or chemical abundances. In particular, time dependence of- 
fers a natural means to sustain strong and broad lines at late times, 
as observed in all Type II SN spectra, without invoking radioactive 
contributions. Such non-thermal excitation/ionization is expected 
to be relevant at the photosphere after a few months, but not as 
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Figure 18. Correction factors versus corresponding color temperature T5 for the 5 1 time-dependent models computed in the present study (dots), and 
shown for bandpass sets BV (left), BVI (middle), and VI (right). For comparison, we also ov erplot crosses for th e {^s^^s} couples obtained with steady- 
state CMFGEN models for SN 1999em jPessart & Hillieill2006j) . SNe 2005cs and 2006bp jPessart et ai]|2007h . There is broad agreement between the 
time-dependent and steady-state predictions, although at low temperatures there is a tendency for the time-dependent correction factors to lie at the lower edge 
of the distribution for the steady-state models. 



soon as a few weeks, and is subject to strong variations accross the 
SN class (SNe 1999br or 2005cs, which boast a strong Hq at the 
recombination epoch, have symptomatically low ^®Ni yields com- 
pared to the standard Type Il-Plateau SN 1999em). 

The time-dependence effects we observe just one week after 
explosion in the outer, fully ionized, SN ejecta also support the idea 
that the energy gain drawn out of recombining ions is not a funda- 
mental driver. Time-dependence effects can lead to over-ionization 
for He, C, N, O, and metals, under fully-ionized hydrogen condi- 
tions. We thus surmise that time dependence should operate, with 
a magnitude to be determined, in SN ejecta of all types, primarily 
because they all combine the properties of fast expansion and low 
density. 

We have investigated the potential impact of time depen- 
dence on the correction factors used in the Expanding Photosphere 
Method, or, essentially, whether the influence on the electron den- 
sity generates a systematic shift in the magnitude of flux dilution. 
We find no sizable and systematic deviation from the correction 
factors obtained, at the same color temperature, with steady-state 
CMFGEN models. This supports our use of steady-state CMFGEN 
models for distance determinations based on early-time Type II SN 
observations (Dessart & Hillier 2006a; Dessart et al. 2007). How- 
ever, the good agreement at late times between time-dependent 
CMFGEN models and late time photospheric-phase observations 
motivates the extension of the time baseline to include the recombi- 
nation epoch, thus allowing for multi-epoch observations that cover 
the entire Plateau phase, thereby reducing the errors on the inferred 
distance. 

For the modeling of the longer evolution of Type II SN ejecta 
and their radiation, a time-dependent approach is warranted. To 
achieve a higher level of consistency, CMFGEN is under develope- 
ment to follow as well the time-dependent evolution of the radiation 
field, together with options for chemical stratification and energy 
deposition from radioactive isotopes. This versatility will allow us 
to study SN ejecta of any type and over months after explosion. 
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Table Al. Term comparison. Note that PR(n=l)/RR is the photoioniza- 
tion/recombination rate from/to level 1 . 
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APPENDIX A: 

We illustrate how the time-dependent terms can influence the He I 
ionization in SN models. We choose a model with a Hubble age of 
6.75 days, and choose a single depth in the atmosphere for which 
helium is still appreciably ionized in the time-dependent model, but 
neutral when such terms are not included. 
In Table |aI] 

we list the various processes, and their rates, 
which influence the helium ionization balance. From this table we 
see the following: The most important process populating level Is^ 
(the ground state) is decays from the 2p ^P° state through the op- 
tically thick resonance transition, followed by the slight imbalance 
between photoionization from level Is"^, and direct recombination 
to level Is^. Many other processes also play a role, including colli- 
sional de-excitation from the 2s & 2p levels, decay from the 2p ^P° 
state, and decays from other members of the np ^P° (n > 2) series. 
As expected, the sum of all these decay processes is the same as 
Dn(He[n=l])/Dt. Note that Dn(He(n=l))/Dt= 1.2 x 10* cm~^ s"^ 
if helium was to change from fully ionized, to fully neutral, over 
one day. Although not directly apparent from the table, most of the 
recombinations to excited states of neutral helium are balanced by 
ionizations from the 2s and 2p states. 
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